Dilute moment n-type ferromagnetic semiconductor Li(Zn,Mn) As 
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We propose to replace Ga in (Ga,Mn)As with Li and Zn as a route to high Curie tempera- 
ture, carrier mediated ferromagnetism in a dilute moment n-type semiconductor. Superior material 
characteristics, rendering Li(Zn,Mn)As a realistic candidate for such a system, include unlimited 
solubility of the isovalent substitutional Mn impurity and carrier concentration controlled indepen- 
dently of Mn doping by adjusting Li-(Zn,Mn) stoichiometry. Our predictions are anchored by detail 
ab initio calculations and comparisons with the familiar and directly related (Ga,Mn)As, by the 
microscopic physical picture we provide for the exchange interaction between Mn local moments 
and electrons in the conduction band, and by analysis of prospects for the controlled growth of 
Li(Zn,Mn)As materials. 

PACS numbers: 75.50.Pp,75.30.Hx,73.61.Ey 



(Ga,Mn)As is a prototypical of a unique class of spin- 
tronic materials in which ferromagnetic coupling between 
dilute local moments is mediated by semiconductor band 
carriersi^ This unusual behavior is realized in the con- 
ventional semiconductor GaAs using only one type of 
dopant, MnQ a , which provides both local spins and holes. 
The simplicity of such a system inevitably brings also 
limitations to the structural, ferromagnetic and semicon- 
ducting properties due to the low solubility of Mn and 
due to the lack of independent control of local moment 
and carrier densities. Among the negative consequences 
are Curie temperatures below room temperature, and p- 
type conduction only. In our study we predict that these 
limitations can be lifted by a straightforward substitution 
of the group-Ill element Ga with group-I Li and group-II 
Zn elements. 

LiZnAs is a stable direct-gap semiconductor which can 
be grown by the high temperature reaction of elemental 
Li, Zn, and As^i Its crystal structure and band struc- 
ture are very similar to those of GaAs. As shown in 
Fig. ^a) and (b), the LiZnAs tetrahedral lattice can be 
viewed as a zincblende ZnAs binary compound, analo- 
gous to GaAs, filled with Li atoms at tetrahedral in- 
terstitial sites near As. The relatively small first ion- 
ization energy of Li makes the (Li) + (ZnAs) _ compound 
half ionic and half covalenti^ The measured band gap of 
LiZnAs (1.61 eV) is very similar to the GaAs band-gap 
(1.52 eV)£ An overall similarity of the LiZnAs and GaAs 
electronic structures, including valence band and conduc- 
tion band dispersions, ground state charge density, and 
phonon dispersion relations has been reported in ab initio 
local-density-approximation (LDA) studiesA^ 

Application of the LDA description of the electronic 
band-structure to Mn-doped zincblende semiconductors 
runs into a conceptual difficulty in dealing with local mo- 
ment levels coincident with itinerant electron bands. To 
partly remedy this inadequacy of standard first princi- 
ples approaches we use in this paper the LDA+U tech- 
nique which combines LDA with the Hubbard descrip- 



tion of strongly correlated localized orbitalsi&2ii2iii The 
LDA+U method is implemented within the framework of 
the first-principles, tight-binding linear muffin-tin orbital 
approach; disorder effects associated with random distri- 
bution of Mn and other defects if present are accounted 
for using the coherent potential approximation (CPA). 12 

In Fig. dc) and (d) we compare LDA+U/CPA den- 
sity of states (DOS) of LiZnAs and GaAs doped with 5% 
of Mn. Since Mn introduces one hole per Mn in GaAs 
but is isovalent in LiZnAs, the structure considered in 
Fig- Did) is Lio.95(Zn .95Mn .o5)As to allow for a direct 
comparison of p-type ferromagnets realized in the two 
hosts. As expected from the similar host band struc- 
tures and nearly atomic-like character of Mn d-states in 
(moderately) narrow-gap zincblende semiconductors, 13 
both systems are ferromagnetic with comparable valence- 
band exchange splittings and similar Mn d and host sp- 
projected DOSs. In particular, the main peak in the 
majority-spin d-orbital DOS is well below the valence 
band edge resulting in states near the Fermi energy hav- 
ing a predominantly As p-character in both systems. 

The LDA+U/CPA band structure of (Ga,Mn)As is 
consistent with a number of experiments in this ferro- 
magnetic semiconductor £ We expect, based on the above 
similarities between the two systems, that Li(Zn,Mn)As 
is also described reliably by this theoretical technique. 
In the following paragraphs we use the ab initio band 
structures to explore Li(Zn,Mn)As mixed crystals over a 
wide range of local moment and carrier dopings, includ- 
ing regimes which are technologically inaccessible in the 
GaAs host. In particular, we focus on the equilibrium sol- 
ubility of isovalent Mn in LiZnAs for concentrations well 
above 1% and on ferromagnetism in LiZnAs with n-type 
conduction which, as we show, can be readily achieved 
without introducing any additional chemical elements. 

The calculated formation energies of Mnz„ are plotted 
in Fig. Ufa) for insulating stoichiometric Li(Zn,Mn)As, 
for p-type Li(Zn,Mn)As with Li vacancies, and for n- 
type Li(Zn,Mn)As with additional Lii atoms occupying 
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FIG. 1: Top panel: Schematics of (Ga,Mn)As and 
Li(Zn,Mn) As crystal structures. Bottom panel: ab initio total 
(blue lines) and Mn d-orbital projected (red filled lines) den- 
sity of states of Gao.95Mno.05As and Lio.9s(Zno.95Mno.o5)As 
mixed crystals. Energy on the x-axis is measured from the 
Fermi energy. 



tetrahedral interstitial sites near Zn. Independent of the 
charge doping, the formation energy of Mnz n is negative, 
i.e., we find no equilibrium solubility limit for substitu- 
tional Mnzn, consistent with its isovalent nature. (Note 
that we focus in this paper on dilute Mn systems far 
from the LiMnAs limit; the LiMnAs compound is an 
antiferromagnelpi^ii^ due to short-range Mn-Mn superex- 
change.) This result contrasts with the positive forma- 
tion energy of MnQ a in GaAsi& which leads to an equi- 
librium solubility limit below 1%. 

A detailed analysis of non-stoichiometric n-type 
Li(Zn,Mn)As is presented in Figs. |^b)-(d). Panel (b) 
demonstrates that in systems with over-stoichiometric Li 
concentrations a large number of Lii atoms can be incor- 
porated at interstitial sites near Zn. Formation energy 
of these single-donor impurities is negative and decreases 
with increasing Mn doping. 

In covalent semiconductors, carrier doping is limited 
due to strong tendency to self-compensation. In the 
LiZnAs host with excess Li, Liz n antisites represent nat- 
ural single-acceptor defects compensating the interstitial 
Lii donors. Ab initio calculations of formation energies 
can be usecUSiii to estimate the dependence of Lii and 
Liz n partial concentrations on the total density of excess 
Li atoms above the Li-(Zn,Mn) stoichiometry. Results 
are plotted in Figs. |2Ic) and (d) for Mn doping of 5% 
and 12.5%, respectively. Although the tendency to self- 
compensation is clearly apparent, the theoretical data 



suggest that large net electron densities are feasible in 
Li(Zn,Mn)As owing to the partly ionic character of the 
compound. Remarkably, the n-type doping efficiency by 
excess Li increases significantly with increasing Mn con- 
centration. 

We now consider the magnetic properties of 
Li(Zn,Mn)As semiconductors. The compatibility of 
the CPA with Weiss mean-field theory allows us to 
estimate the strength of Mn-Mn magnetic coupling, 
at a given chemical composition, from the calculated 
energy cost, E rev , of reversing one Mn moment with 
all other Mn moments alignediiS*^ Results for E rev 
shown in Fig. |3^a) were obtained using a rigid-band 
approximation. Here the LDA+U/CPA band structure 
was calculated in one charge state of the compound only, 
and the dependence on carrier doping was obtained by 
shifting the Fermi energy while keeping the bands fixed. 
Note that for the 5% Mn doping considered in the figure, 
carrier doping of —5% corresponds to one hole per Mn 
and +5% to one electron per Mn. Open and closed 
squares represent E rev obtained applying the above rigid 
band scheme, based on n-type Lii.o5(Zn .95,Mn .o5)As 
and p-type Lio.95(Zno.95,Mno.o5)As band structures, 
respectively. The similarity of the two curves indi- 
cates that the position of the Fermi energy plays the 
dominant role in magnetic interactions while disorder 
effects associated with non-stoichiometric Li-(Zn,Mn) 
configurations are less important. For comparison we 
plot in Fig. |HIa) the rigid band E rev obtained from 
the Gap.95Mno.05As ab initio band structure (see also 
Ref. |20|) . As expected the two compounds show very 
similar magnetic behavior. 
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FIG. 2: Ab initio impurity formation energies of Li(Zn,Mn) As 
mixed crystals: (a) Mnzn formation energy as a function of 
Mnzn doping for stoichiometric structures (red open dots) , for 
p-type systems with 5% of Li vacancy acceptors (black trian- 
gles), and for n-type systems with 5% of donor Lii impurities 
occupying Zn-tetrahedral interstitial sites (blue diamonds), 
(b) Lii formation energy near the stoichiometric composition 
as a function of Mnzn doping, (c) Partial concentration of Lii 
(blue diamonds) donors, Lizn acceptors (red dots), and net 
electron doping (black triangles) for 5% of Mnz n as a func- 
tion of total concentration of excess Li (Lii plus Liz n ). (d) 
Same as (c) for 12.5% of Mn Zn . 
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FIG. 3: Energy cost of nipping one Mn moment calculated 
with all other Mn moments aligned. Open and closed squares 
were obtained from n-type Lii.05 (Zno.95,Mno.os)As and p- 
type Lio.95(Zno.95,Mno.o5)As ab initio band structures, re- 
spectively (see text for details). Circles represent calculations 
for Gao.95Mno.05As. 




5 10 12 3 



Carrier density (10 20 cm" 3 ) d/a 

FIG. 4: (a) Model calculations of Mn&i d-orbital hybridiza- 
tion potentials with GaAs sp-orbitals at the Fermi energy in 
the valence band (blue diamonds and green triangles) and 
conduction band (red dots) as a function of hole density and 
electron density, respectively, (b) Model band structure of 
the GaAs host illustrating the filled conduction band states 
(red full line) and empty valence band states ( black and blue 
lines) for electron or hole density of 10 21 cm" 3 (largest car- 
rier density considered in panel (a)), (c) Ab initio Mn-Mn ex- 
change interactions as a function of Mn interatomic distance 
evaluated along the [110] crystal direction. 



The calculated positive values of E rev at large hole 
concentrations correspond to ferromagnetic Mn-Mn cou- 
pling, which is consistent with experimentally observed 
ferromagnetic ground states in highly doped p-type 
(Ga,Mn)As. (At low enough hole concentrations, E rev 
changes sign and fcrromagnetism becomes unstable due 
to the dominant role of short-range antifcrromagnetic 
Mn-Mn superexchangei^) The comparable strength of 
ferromagnetic Mn-Mn coupling for n-type and p-type 
doping at high carrier concentrations suggests that ex- 
change interactions between Mn (i-states and conduction 
band states in zincblende semiconductors can be strong 



and that Li(Zn,Mn)As with properly adjusted stoichiom- 
etry is a high Curie temperature n-type ferromagnetic 
semiconductor. 

We now discuss the physical origin of the high temper- 
ature electron-mediated fcrromagnetism. In zincblende 
semiconductors doped with Mn d 5 impurities, the local 
moments can have either a direct exchange interaction 
with band electrons on the same site and/or an interac- 
tion due to sp—d hybridization between the Mn local mo- 
ment and band electrons on neighboring sitesi2i*22i The 
latter interaction is much stronger and high-temperature 
ferromagnetism is expected to occur only when the hy- 
bridization coupling mechanism is present. In GaAs, the 
valence band edge at the T-point has a p-orbital charac- 
ter with a larger contribution from As^ Since the den- 
sity of states at the valence band edge is large, the oc- 
cupied states for typical hole densities in ferromagnetic 
(Ga,Mn)As are still close to the T-point and have sim- 
ilar atomic orbital composition. The p — d hybridiza- 
tion is allowed by symmetry within the whole Fermi sea 
(including the T-point) which, together with the large 
magnetic susceptibility of heavy mass holes, explains the 
experimentally observed robust ferromagnetism in p-type 
(Ga,Mn)As. 

The bottom of the conduction band of GaAs has an s- 
orbital character with a larger contribution from Ga. The 
smaller admixture of orbitals from the nearest neighbor 
As atoms to MnQa, the vanishing s — d hybridization at 
the T-point due to symmetry, and the small magnetic sus- 
ceptibility of low effective mass electrons near the bottom 
of the conduction band all suggest that ferromagnetism 
is unfavorable in n-type systems. As shown in the inset 
of Fig. |21 a weaker tendency towards ferromagnetism in 
n-type materials is indeed observed in our ab initio cal- 
culations, however, is limited to relatively low electron 
dopings. 

To explain the seemingly unexpected robust ferromag- 
netism at larger electron concentrations we present, in 
Fig-EI a ) and (b), tight-binding model calculations of the 
energy dependent sp — d hybridization potential in the 
conduction band and compare with the more familiar 
valence band case. When moving off the T-point, the 
hybridization rapidly sets in partly due to the s — d hy- 
bridization allowed by symmetry at non-zero wavevectors 
and partly due to the admixture of As and Ga p-orbitals 
which increases with increasing wavevector. The strength 
of the sp — d hybridization in highly doped n-type sys- 
tems becomes comparable to that of p-type materials as 
the Fermi energy leaves the bottom part of the conduc- 
tion band and approaches the L-point. The same picture, 
which we have discussed for the more familiar GaAs host 
doped with Mn, applies to Li(Zn,Mn)As whose material 
properties, as analyzed above, are much more favorable 
for achieving large n-type doping. 

Within the mean-field theory, the ferromagnetic Curie 
temperature, T c , can be estimated as T c ss E rev /&kB, 
where ks is the Boltzmann constant. The mean-field 
approximation should be reliable when the carrier medi- 
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ated Mn-Mn coupling is sufficiently long range but tends 
to overestimate T c when the carriers become more lo- 
calized and magnetic interactions short-ranged £12&&L In 
Fig. 0Jc) we plot the spatially dependent interatomic ex- 
change energies obtained by mapping the LDA+U/CPA 
total energies on the Heisenberg Hamiltonianpi^ As seen 
in the figure, the n-type Li(Zn,Mn)As shows very sim- 
ilar magnetic interaction characteristics to the p-type 
(Ga,Mn)As with corresponding doping densities. In both 
systems the leading exchange interactions are ferromag- 
netic and the interaction range safely exceeds the average 
Mn-Mn moment distance. Consistent with the theory ex- 
pectations, experimental high quality (Ga,Mn)As mate- 
rials show mean field-like magnetization curves and Curie 
temperatures proportional to the density of Mnc a ^ 
Maximum Mnc a doping achieved so far is approximately 
6% and the corresponding record T c is 173 Kii£ Calcula- 
tions shown in Fig. 0] suggest comparable Curie temper- 
atures for the n-type Li(Zn,Mn)As counterparts. Since 
we found no doping limit for Mnz n in Li(Zn,Mn)As and 
straightforward means of high electron doping in the ma- 
terial, Li(Zn,Mn)As might lead not only to the realiza- 
tion of a high Curie temperature diluted magnetic semi- 
conductor with electron conduction but also might allow 
for larger magnetic moment density and therefore larger 
maximum T c than (Ga,Mn)As. 

We conclude by briefly discussing prospects for the 
controlled growth of Li(Zn,Mn)As materials. Bulk high 
structural quality LiZnAs crystallites have been formed 



by reaction of near equimolar elemental Li, Zn, and As>^i 
The crystallites are stable in inert atmospheres but slowly 
oxidize in air. It is probable that Li(Zn,Mn)As crystal- 
lites could be formed in the same way. However, to obtain 
the required control of stoichiometry to produce the fer- 
romagnetic semiconductors predicted by our calculations, 
molecular beam epitaxy (MBE) is likely required. LiZ- 
nAs has not previously been grown by MBE, presumably 
because non-magnetic LiZnAs does not seem to have any 
advantage over GaAs. The similarity of Li(Zn,Mn)As to 
(Ga,Mn)As suggests that it should be possible to grow 
Li(Zn,Mn)As by MBE; here Li can be supplied from a 
metal-organic gas source. Li(Zn,Mn)As could be grown 
by MBE on lattice matched, relaxed, (In,Ga)As epilayers 
on normal GaAs substrates, carrier doping achieved by 
careful control of the Li flux, and, furthermore, it should 
be possible to avoid post growth oxidation by overgrowth 
of the Li(Zn,Mn)As epilayers. 
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